First observation of the decay K + 
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Experiment 865 at the Brookhaven AGS has observed the decay K + — » e + v e n + . The branching 
ratio extracted is (1.72±0.37(stat)±0.17(syst)±0.19(model)) x 10 -8 where the third term in the error 
results from the use of a model to extrapolate into a kinematic region dominated by background. 

PACS numbers: 13.20.Eb, 13.40.Ks 



> 



o 
in 
o 
in 
o 

">< 

<D 
Oh! 



X 



The internally converted, radiative K%2 decays, K + — » 
l+vl' + l'~ , are an important source of information on the 
kaon. For example, within the framework of Chiral Per- 
turbation Theory (ChPT) |l( radiative kaon decays can 
serve both as an important test and a source of input 
parameters for the theory. 
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FIG. 1: Graphs for contributions to K + — > e + v/j, + /i~ . 

The K + — » l + vl l+ V is described by the graphs of 
Fig. ^ The tree diagrams (a) and (b) are the inner 
bremsstrahlung (IB), where the virtual photon is ra- 
diated by the kaon or positron. This contribution is 
electron-helicity suppressed and is negligible for the de- 
cay K + — > e + vfi + ii~. The short distance, structure- 
dependent (SD) terms are combined in graph (c). The 
SD contribution is characterised by form factors Fy, Fa, 
and R, which we define in accordance with the Particle 
Data Group ,2j . These may be functions of W 2 and q 2 , 
the invariant masses of the e + v and the pairs, re- 

spectively. In the vector meson dominance picture 
this dependence has the form 



branching ratio of 1.12 x 10 8 . The previous experimen- 
tal limit was < 5.0 x 10" 7 0,0. 

Experiment E865 at the Brookhaven National Lab- 
oratory Alternating Gradient Synchrotron (AGS) has 
produced substantial improvements in our knowledge of 
these radiative decays. Results for the decays K + — > 
H + ve + e~ and K + — > e + i/e + e~ have already been re- 
ported [q . This paper presents results from E865 for the 
first observation of the decay K + — > e + v/i + /i~ . 
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FIG. 2: Plan view of the experimental apparatus for E865. 



Fy(q 2 , W 2 ) = F v (0, 0)/[(l - q 2 /m 2 p )(l - W 2 /m 2 v )} (1) 

with similar expressions for Fa and R. Here, m p is the 
p-meson mass and my is the mass of the if* (892) for Fy 
and of the if* (1270) for F A and R. ChPT relates the 
form factors Fy and Fa to those of the it —* evj decay 
and form factor R to the kaon charge radius. Bijnens el 
al. 0] gave a ChPT prediction for the K + — » e + vji + fj," 



The experimental apparatus for Brookhaven E865 is 
shown in Fig. [2 and has been described in detail else- 
where [jj. A 6 GeV/c unseparated beam from the AGS 
was incident on a 5-m long evacuated decay volume. De- 
cay particles from this region were separated into posi- 
tive and negative charges by a dipole magnet and were 
momentum-analysed by a spectrometer containing a sec- 
ond dipole magnet and four 4-view wire chambers, PI 
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- P4. Particle identification was provided by two pairs 
of gas Cherenkov counters, CI and C2 (hydrogen on 
the negative-particle side and methane on the positive 
side), a 30 x 20-element Shashlyk calorimeter containing 
15 radiation lengths of lead/scintillator sandwich, and 
a muon detector containing twelve layers of iron with 
twelve 2- view wire chambers interspersed. Additionally, 
there were 4 hodoscope planes, A - D, for timing and 
triggering. A 12 (horizontal) x 32(vertical)-element pixel 
detector, with pixel size 7x7 mm, was located in the 
incident kaon beam to determine the position of the de- 
caying kaon at the entrance to the decay volume. 

The data taking for this part of E865 took place in 
parallel with studies of the K + — > ir + n~e + v (K^a) and 
K + —>■ 7t + u + m~ decays, which have already been pub- 
lished US, 03. 

Selection of candidate events required three tracks giv- 
ing a vertex z-coordinate within the decay region and an 
acceptable value of S, where S 2 is the sum of squares 
of the deviations of the three tracks from the fitted ver- 
tex. Also, the individual tracks were each required to 
have good \ 2 f° r the reconstruction and good timing. 
Electron identification required signals in both positive- 
side Cherenkov detectors and an energy deposition in the 
calorimeter equal to the reconstructed track momentum. 
The above cuts were also used in the event selection for 
the K e 4 analysis. Additionally, for K + — » e + z^/i + /i~, the 
muon candidate tracks were required to have sufficient 
hits in the muon wire chambers, a hit in the appropri- 
ate element of the B hodoscope, located in the middle of 
the muon stack, and a signal consistent with minimum 
ionising in the shower calorimeter. 

For each event, the neutrino momentum was calcu- 
lated from the missing momentum; p„ = Pk+ — — 
p M — p e + . The magnitude of Pk+ was taken as the aver- 
age value, derived from measurements of pion momenta 
from decays. The direction of Pk+ was determined 
from the beam pixel detector and the reconstructed ver- 
tex where possible. About 45% of events had an unam- 
biguous hit in the beam pixel detector. For those events 
that did not, the average kaon beam direction was as- 
sumed in the event reconstruction. 

Simulation of the experiment was carried out using the 
GEANT package A problem with this package is 

that pion interactions are not always well simulated. As a 
result, the probability that a pion can penetrate well into 
the muon stack, causing it to be misidentified as a muon, 
is not well determined by the simulation. Experimental 
studies using K + — > 7r + 7r°, K + — > ir n + v (each followed 
by 7T° — > 7e + e~) and K + — > 7r + 7r + 7r~ events established 
that this 7T to \x misidentification probability is calculated 
by the simulation with an uncertainty of about 10%. 

With the above cuts, 1834 candidate events remain. 
These are shown in Fig. as a scatter plot of m e + M + M - 
against m e + V n+u--, which also shows the region of 
TO e +j, M + M - where genuine events are expected. 
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FIG. 3: Scatter plot of m e + li + IJi - against m e + yM + M - for can- 
didate K + —* e + ffi + fi~ events. Valid events lie in the band 
m e+t//j+n- ~ mK indicated on the plot. 

Three sources of background were considered: 

(a) K + — > ix + ii~e + v (K e 4) with both pions 

misidentified as muons, 

(b) Accidentals, and 

(c) K + — » 7r + 7r + 7r _ (A'3^) with the pions 

misidentified as e + , fi + and fjT respectively. 
The contributions from these were determined as follows. 

K e 4 events give the largest contribution to the back- 
ground. A large sample of Monte Carlo events (6 x 10 7 , 
equivalent to 20 times the number of data events) was 
generated. Background hits taken from actual data 
events were added to the simulated events, and the sam- 
ple was analysed as for data events. This procedure 
should give a realistic estimate of the background shape. 
The magnitude of he K e 4 background has an error due 
to uncertainty in the tt to [i misidentification probability, 
which is about 0.05 with an error of ~ 10%. 

A sample of accidental events was extracted from the 
data by selecting events with bad timing and badly re- 
constructed vertices. The reconstructed total charged- 
track momenta for these events shows a tail above 6.5 
GeV/c. The total charged-track momentum spectrum 
for the K + — > e + v^i + ^ candidates shows a similar 
tail, which is assumed to arise entirely from accidentals. 
Therefore, the background contribution from accidentals 
was assumed to have the same shape as the bad-timing, 
bad-vertex events, with a magnitude derived by scaling 
this high-momentum tail to match that in the data. 

The background contribution from decays was es- 
timated by examining the K + — » 7r + 7r + 7r~ data events 
obtained from the minimum-bias trigger. The fi + and /i~ 
identification was required but there was no particle iden- 
tification requirement on the third particle. The sample 
was then scaled by the known tt + to e + misidentification 
probability of 1.3 x 10~ 3 , determined from a study of 
A e 4 events, where the contamination from decays is 
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FIG. 4: Simulation results for the K + — > e + vfi + pT (top left- 
hand plot) and the computed contributions from the three 
background sources, K e 4, and accidentals. 



easily identified by the kinematics of the 3ir final state. 

Background from K + — > 7r + 7r _ /x + f is negligible be- 
cause of the small 7r + to e + misidentification probability 
(1.3 x 10~ 3 ) and the low branching ratio for this decay. 

The three sources of background are shown as plots of 
m e+fj,+fi- against m e + u ^+^- in Fig. 0] which also shows 
the simulated K + — > e + vp + fi~ signal. The signal ap- 
pears as a peak in m el/flfl at the kaon mass, 0.494 GeV/c 2 , 
the region indicated by the vertical lines in Figs. and 
0] It is also apparent from Fig. that the signal-to- 
background ratio will be best for high values of m eW j. 

Because the signal-to-background ratio becomes poor 



at low m eMM , it is difficult to extract significant informa- 
tion on form factors from the present experiment. Ini- 
tially, therefore, we assumed form factors for K + — > 
e + vfi + pT in order to extrapolate the signal into regions 
where it is not observable, to extract the total branching 
ratio. To do so, the 2-dimensional plot, Fig. EI was fitted 
with the sum of the signal and background contributions 
shown in Fig. 0] While the magnitude of the signal was 
a free parameter in the fit, the shape of the signal dis- 
tribution was fixed by the form factors Fa, Fy and i?, 
whose values are taken from the E865 measurements on 
K + — > e + ve + e~ and ii + ve + e~ 0: 

PV(0,0) = 0.112 ± 0.015(stat) ± 0.010(syst) 

±0.003(modcl), (2) 



Fa (0,0) = 0.035 ± 0.014(stat) ± 0.013(syst) 

±0.003(model), (3) 
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FIG. 5: Projection onto the m el , w axis of the data of Fig. H 
(points) and of the fitted distribution (unshaded histogram). 
The three background contributions are shown separately 
(shaded histograms). 



R(0,0) = 0.227 ± 0.013(stat) ±0.010(syst) 

±0.009(model). (4) 

Since the K e ± contribution to the background is uncer- 
tain to about 10%, we allow the magnitude of it to float 
in the fitting. 

The fit is shown in Fig.[3]as a projection onto the m eiw 
axis of the fit and of the data from Fig. E] The three 
contributions to the background are shown separately. 
Fig.Elshows the same projection for events in three bands 
of 

^e/x/i as follows: 

Fig. EJa): 0.2825 < m ew < 0.355 GeV/c 2 , 

Fig.Etb): 0.355 < m eWi < 0.4275 GeV/c 2 , 

Fig. Etc): m ew > 0.4275 GeV/c 2 . 
These regions are indicated by the horizontal lines in 
Figs.Eland^] Fig. Etc) contains 24 data events with a fit- 
ted background of 6.9 events, and Fig.Etb) has 176 events 
with a background of 133±4 events for the reconstructed 
kaon mass range 0.465 < m eJW < 0.540 GeV/c 2 . 

The fit resulted in a branching ratio, B, of 

B = (1.72 ±0.37(stat)±0.17(form factor) 

± 0.09(slope) ± 0.17(syst)) x 10" 8 (5) 

The K e 4 decay was used as a normalisation channel for 
the analysis. Thus the systematic error arises predomi- 
nantly from uncertainties in simulation of the muon de- 
tection efficiency since muons are not involved in K e ^. 
For this, we estimate ±10%. Additional errors arise from 
uncertainties in the form factors of Eqs. J5J, J2J and 
and also in the slopes of the q 2 and W 2 dependence of 
these form factors. For this contribution, we assume, as 
m Ref. [|, that these slopes are uncertain to ±30%. The 
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FIG. 6: Distribution 



of events in three bands of m e 



0.355 GeV/c 2 < 



< 0.4275 GeV/c 2 and (c) 



(a) 0.2825 GeV/c 2 < 




evun (GeV/c^ 

0.355 GeV/c 2 , 



(b) 



> 0.4275 GeV/c 2 . The points are the experimental data and the 



histogram is the fitted spectrum. The shaded areas show the fitted background contribution. 



fit gave x 2 — 132 for 112 degrees of freedom. The best 
fit resulted in a scaling factor for the K e \ background of 
0.86 ± 0.02 which is consistent with the estimated uncer- 
tainty in the 7r to fi misidentification probability of about 
10%. Combining the errors in quadrature gives 

B = (1.72 ± 0.45) x 10~ 8 . (6) 

A study of the statistical significance of the background 
function shows that, in fitting the entire m ejlil vs. m evfJjfll 
distribution, the probability of a statistical fluctuation in 
the background simulating the signal is < 10~ 6 . 

If Fa, Fy and R are all allowed to vary in the fitting, 
the resulting values for B, Fa, Fy and R are consistent 
with Eqs. ©, ©, ® and © but with substantially 
larger errors. This is the expected consequence of the 
variation of signal-to-background ratio across the range 
of phase space covered by our data. However, since R is 
the dominant term, a fit was carried out with R varied, 
with Fa and Fy constrained at the values of Eqs. J2J) and 
©. This fit gave 

R(0, 0) = 0.303 ± 0.036(stat) ± 0.011 (form factor) 

± 0.009(slope) ± 0.017(syst). (7) 

This is consistent with Eq. J3Jl at the 1.7-standard- 
deviation level, and although our results cannot im- 
prove our knowledge of the form factors, they can at 
least demonstrate consistency with the values from the 
K + — > l + ve + e~ data. 

The result for B can be compared with the ChPT pre- 
diction of Bijncns et al. 4] of 1.12 x 10~ 8 . If the form 
factors of Eqs. J5J), @ and J3J are used in the theo- 
retical calculation, the prediction becomes B — (1.04 ± 
0.15) x 10~ 8 which is consistent with our result at the 
1.5 standard-deviation level. 

In summary, we have made a first observation of the de- 
cay K + — » e + is[i + fi~ and have determined the branching 
ratio to 25%. We find values for the form factors that are 



consistent with those for other K + — > l + vl l+ V decays 
but with less accuracy. The branching ratio is reasonably 
consistent with a ChPT prediction. 
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